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Abstract. One of the most intriguing discoveries made by the Swift satellite is the flaring activity in about half of 
the afterglow lightcurves. Flares have been observed on both long and short duration GRBs and on time scales 
that range from minutes to ~1 day after the prompt emission. The rapid evolution of some flares led to the 
suggestion that they are caused by late central engine activity. Here, I propose an alternative explanation that 
does not need reviving of the central engine. Flares can be powered by delayed magnetic dissipation in strongly 
magnetized (i.e. with initial Poynting to kinetic flux ratio >1) ejecta during its deceleration due to interaction with 
the external medium. A closer look at the length scales of the dissipation regions shows that magnetic dissipation 
can give rise to fast evolving and energetic flares. Multiple flares are also expected in the context of the model. 
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1. Introduction 

A rich early afterglow phenomenology has recently been 
revealed thanks, to a large extent, to the X-ray detector 
on board to the Swift satellite. The systematic study of X- 
ray afterglows (Zhang et al. 2006; Nousek et al. 2006) has 
shown that they are typically characterized by an initial 
steep decay, associated with the end of the prompt emis- 
sion phase, followed by a flatter decay that lasts for about 
one hour and probably corresponds to the onset of the af- 
terglow. At later times, the lightcurves are characterized 
by the "normal" decay observed before Swift. 

Another exciting Swift discovery is the afterglow flares 
observed in about half of the bursts (Burrows et al. 2005). 
Flares appear on different timescales if after the prompt 
emission ranging from ~ 100 to ~ 10 5 sec, they are char- 
acterized by duration Stf such that 5t{/t{ < 1, and show 
spectral evolution with respect to the smooth decay part 
of the afterglow. The flares are involving substantial en- 
ergy release that ranges from a small fraction up to an 
amount equal to that of the prompt GRB emission (as in 
the case of GRB050502b; Falcone et al. 2006) and they 
have been observed in the afterglow lightcurves of both 
long and short duration GRBs (Campana et al. 2006). 

Although bumps in the afterglow lightcurves can be a 
result of slower shells catching up with the decelerating 
ejecta (e.g. Rees & Meszaros 1998) or interaction of the 
ejecta with a clumpy external medium, the rapid evolution 
of some of the flares has been used as an argument against 
their external shock origin. This led to the suggestion that 
the central engine revives at later times (ranging from 
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hundreds of seconds to days) giving rise to flares through 
late internal shocks or magnetic dissipation (Burrows et 
al. 2005; Zhang et al. 2006). Suggestions on how to revive 
the central engine have been made by King et al. (2005), 
Dai et al. (2006), Fan et al. (2006), Perna et al. (2006). 

Here, I make an alternative suggestion which does not 
need late time activity of the central engine. I investigate 
the possibility that the deceleration of strongly magne- 
tized ejecta leads to revival of MHD instabilities, that were 
delayed by time dilation before the deceleration phase. 
The resulting dissipation produces flares in localized re- 
connection regions in the flow. 

2. Magnetic dissipation in a decelerating shell 

Magnetic fields may be the main agent extracting en- 
ergy from the central engine, leading to a Poynting-flux- 
dominated flow (e.g., Thompson 1994; Meszaros & Rees 
1997; Spruit et al. 2001). Magnetic dissipation through in- 
stabilities in an axisymmetric flow (Lyutikov & Blandford 
2003; Giannios & Spruit 2006) or directly through recon- 
nection in a highly non-axisymmetric flow (Drenkhahn & 
Spruit 2002) can power the prompt emission with high 
radiative efficiencies and accelerate the flow to ultrarela- 
tivistic speeds. The photospheric emission of a dissipative 
MHD flow is highly non-thermal and can be responsible 
for the prompt GRB emission (Thompson 1994; Giannios 
2006). 

Dissipation of Poynting flux is partial, however, if it de- 
pends on global magnetic instabilities. MHD instabilities 
(such as the kink instability) grow on the Alfven crossing 
time in the fluid frame. Time dilation in the central engine 
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frame by the bulk Lorentz factor of the flow T results in 
slowing down of the dissipation in ultrarelativistic flows. 
Giannios & Spruit (2006) have shown that for large initial 
Poynting to kinetic flux ratios, the flow remains magnet- 
ically dominated up to radii >10 16 cm, i.e. close to the 
distance where the deceleration of the flow against the 
external medium is expected. 

The deceleration of the flow can, however, naturally 
lead to reviving of the instabilities in the regions which 
come into causal contact again and to substantial mag- 
netic dissipation in the afterglow phase. Another (and 
probably complementary) possibility is that dissipation of 
magnetic energy is triggered by the crossing of the reverse 
shock through the magnetized shell (Thompson 2005). 

Here, I investigate the possibility that late time dis- 
sipation takes place in this way and explore whether it 
can power the commonly observed flares in the GRB af- 
terglows. I discuss in particular the energetics and the 
variability expected from such events. 

2.1. Deceleration of a magnetized shell 

The issue of flow acceleration and of prompt emission in a 
Poynting flux dominated flow has been studied in a num- 
ber of works in the context of GRB flows (e.g., Drenkhahn 
& Spruit 2002; Vlahakis & Konigl 2003; Giannios & Spruit 
2006). At a sufficient distance from the source, the flow has 
accelerated to its terminal Lorentz factor and the magnetic 
field is dominated by its toroidal component. After the ul- 
trarelativistic flow has swept up enough mass from the ex- 
ternal medium (which may be the interstellar medium or 
the wind material), it starts decelerating, driving a shock 
into the external medium. In this work, I focus on this 
deceleration phase. 

I assume that a cold, strongly magnetized shell with 
a dominant toroidal magnetic field is moving with initial 
bulk Lorentz factor Po S> 1 into the external medium. 
The magnetic content of the ejecta can be conveniently 
parameterized by the magnetization a defined as the ratio 
of Poynting to kinetic energy flux a = B' 2 /4irw, where 
B' and w = pc 2 + e + p are the magnetic field and the 
enthalpy density of the ejecta respectively as measured by 
an observer comoving with the flow. The enthalpy density 
consists of the rest mass energy density pc 2 , the energy 
density e and the gas pressure p. Since the shell is assumed 
to be cold initially, its magnetization is ao = B' 2 /Anpc 2 . 

The initial phase of the interaction of the magnetized 
shell with the external medium has been studied by Zhang 
& Kobayashi (2005) who have shown that, for "typical" 
GRB parameters, a reverse shock forms in the ejecta as 
long as do 100. In the shocked part of the ejecta, both 
gas and magnetic field contribute to the total pressure; 
their sum balances the pressure of the shocked external 
medium at the contact discontinuity. For high enough ini- 
tial magnetization of the ejecta ao > 0.1, the magnetic 
pressure dominates the gas pressure in the shocked ejecta 
irrespective of whether the reverse shock is Newtonian, 



mildly relativistic or relativistic. This can be shown by 
solving for the MHD shock conditions (Kennel & Coroniti 
1984; Zhang & Kobayashi 2005), where it is also shown 
that shock compression leads to an increase of the mag- 
netization of the ejecta, i.e. cr s h <; <Jq. 

Thus, for a > 0.1, one can neglect the gas pressure of 
the shocked ejecta with respect to the magnetic pressure. 
The pressure balance at the contact discontinuity yields 
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where B'^, r s h are magnetic field strength and the bulk 
Lorentz factor of the shocked ejecta and p c is the density 
of the external medium. From this point on, I focus on the 
shocked ejecta and omit the subscript "sh" for simplicity 
in the notation. 

2.2. The magnetic dissipation region 

I have argued that the deceleration of the flow (and/or 
the crossing of the reverse shock) revives MHD instabil- 
ities that lead to dissipation of magnetic energy through 
reconnection of magnetic field lines at different locations 
in the flow. The energy released in a dissipative event and 
its observed duration depend on the characteristic dimen- 
sions of the dissipation region. Let l[ and 1' 2 be the char- 
acteristic length scales in the radial and perpendicular to 
the bulk motion directions (measured by a comoving ob- 
server) of this region respectively. Furthermore, I consider 
l'i to be the scale over which the field changes polarity 
and that magnetic field lines are advected in the middle 
plane of this region where they reconnect. The reconnec- 
tion plane is assumed to be perpendicular to the radial di- 
rection. More generally, one should consider an arbitrary 
angle between the radial direction and the normal to the 
reconnection plane. However, I have checked that such a 
generalization only complicates the analysis that follows, 
without introducing any essential changes on the results, 
and has been avoided. 

The electromagnetic energy that is contained in the 
dissipation region is (in the central engine frame) 
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where eq. Q has been used in the last step. If a large 
fraction of the dissipated energy leads to fast moving elec- 
trons, they can efficiently radiate most of it if the elec- 
trons' cooling timescale is shorter than the expansion time 
scale £ C xp = r/Tc of the flow. Here, I show that this is the 
case in the dissipation region under consideration. 

An efficient cooling mechanism in the strongly magne- 
tized flow is synchrotron cooling with a time scale 



t syn = 3m c c/4cr T /3g7eZ7 B 



(3) 



where 7 e is the characteristic Lorentz factor at which 
electrons are accelerated in the reconnection region and 
Ub — B' 2 /8ir is the magnetic energy density. The avail- 
able magnetic energy per proton (and electron) is >m p c 2 
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for a flow with <jq <; 1, which is the range of ao that I 
focus on in this study. Dissipation of most of this energy 
to the electrons leads to 7 C > m p /m c . By using, for an 
order of magnitude estimate, T — 50, r = 3 ■ 10 17 cm, 
p c = 1CT 24 gr cm~ 3 , and j c = 10 3 5 in eqs. JIJ and J3J), 
I find t syn /t exp ~ 10 -2 and the characteristic synchrotron 
photon energy E syn — 3Tj 2 q c hB' /2m e c ~ 100 eV, i.e. in 
the soft X-rays. Smaller values for t syn /t cxp and higher for 
-Esyn are found for "typical" parameters that correspond 
to the case where the external medium is the stellar wind. 

Thus, it is possible that a large fraction of the recon- 
nected energy is promptly radiated away in the X-rays 
through synchrotron emission. This energy is emitted in 
the 9 ~ 1/r forward cone, i.e., with a beaming factor 
47r/27r(l — cos#) ~ AT 2 . The beaming leads to isotropic 
equivalent energy Ef° that is AT 2 times larger than the 
dissipated energy 

Ef° = AT 2 E A . (4) 

If the released magnetic energy is to be observed as 
a flare, the dissipation must take place sufficiently fast. 
This leads to conditions on the length scales l[ and l 2 
that have to hold. These conditions are investigated in 
the next section. 

2.3. Variability constraints 

In the dissipation region, magnetic fields are advected 
along a distance ~ l[ into the reconnection plane which 
has a characteristic length 1' 2 . The observed duration of 
the flare Stf is the sum of two contributions. The first is 
the "radial" duration St T coming from the fact that while 
dissipation takes place, the flow expands radially resulting 
at different arrival times for photons that are emitted at 
different stages of the dissipation event. The second con- 
tribution comes from the "angular" duration St ang that is 
related to the delay in arrival of photons that are emitted 
in different locations on the reconnection plane. The du- 
ration of the flare is Stf = St r + St ang . I turn to these two 
contributions and study them separately. 

The "radial" duration of the flare depends on how fast 
reconnection takes place. The speed v r at which magnetic 
reconnection proceeds in a strongly magnetized plasma 
is enhanced by the relativistic kinematics of the problem 
(e.g., Lyutikov & Uzdensky 2003; Lyubarsky 2005) and is 
a substantial fraction e ~ 0.1 of the Alfven speed which 
is close to the speed of light for <7o ^ 1. Magnetic field 
lines are advected for a distance l[/2 above and below 
the reconnection plane with speed w r , which leads to a 
comoving dissipation timescale t' d ~ l' 1 /(2v I ) ~ l[/(2ec). 
Dividing this timescale with the expansion timescale r/Tc, 
one has an estimate on the fractional duration St T /tf of 
the flare because of the radial motion of the reconnecting 
region 

5t r /t { = n[/2er. (5) 

The "angular" duration is related to the characteris- 
tic length V 2 of the reconnecting region. A first constraint 
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to l 2 comes from causality arguments that yield l 2 < r/T. 
The observer sees an emitting region with length V 2 from 
a characteristic angle ~ 1/r. The difference in the arrival 
time of two photons coming from the two edges of the 
emitting region is St ang ~ l' 2 /cT. How the observer time is 
related to the radius and bulk Lorentz factor of the blast 
wave depends on the details of the deceleration profile 
which, in turn, depend on the density profile of the exter- 
nal medium and on the magnetic content of the ejecta. To 
keep the study general, I assume that tf = r/acT 2 , where 
a is a parameter with typical values ~ a few [e.g., for bary- 
onic ejecta decelerated in constant density medium a = 4 
(Waxman 1997) and in stellar wind a — 2 (Pe'er & Wijers 
2005)]. Using the expressions for 5t ang and tf, I have 

St ang /t f = aTl' 2 /r. (6) 

The observed duration of the flare 5tf is given by the 
sum of the radial and angular timescales. Using eqs. JSJ 
and ©, I arrive at 

8tt/tt = Tl'j2er + aTl' 2 /r. (7) 

The observed fractional duration of the flare thus con- 
strains both characteristic length scales of the reconnec- 
tion region 

l[ < 2er(6tf/tf)/T, (8) 
l 2 <r(Stf/tf)/aT. (9) 

In the process of deriving the angular time scale, I have 
implicitly ignored fast fluid motions in the frame moving 
with the bulk of the flow. The material leaves the recon- 
nection region with the Alven speed (e.g. Lyubarsky 2005) 
however, which is close to the speed of light. This can af- 
fect the estimate (JHJ) by allowing for larger values of 1' 2 
(and ultimately for larger energy release in a flare of a 
specific Stf/tf, see the rest of the section). This effect is 
particularly pronounced in uq 1 flows which can con- 
tain emitting regions with high "internal" Lorentz factors 
and lead to more rapid variability than that estimated 
here (Lyutikov & Blandford 2003; Lyutikov 2006). 

Expressions (JHJ and ©, combined with eqs. (@J and 
@, yield a limit on the isotropic equivalent energy for a 
flare that can be produced by a single reconnection event 

Ef° < 20e(5tf/tf) 3 r 2 r 3 p c c 2 /a 2 , (10) 

where the numerical factor of the last expression 64/3 was 
set to 20. Eq. 1)10(1 can be rewritten in a more compact 
form by using the fact that 47rr 3 p e /3 — M swopt , where 
M swcp t is the mass of the external medium that has been 
swept up by the blast wave. Furthermore, the forward 
shock conditions dictate that r 2 M swept c 2 is to be identi- 
fied with the energy that has been passed onto the forward 
shock Elf*. Expression (|10J) can, thus, be rewritten 

£^ so <5e(^ f /i f ) 3 ^ s so /a 2 . (11) 

From the last expression, it is clear that the model pre- 
dicts that the energy available to power the fastest evolv- 
ing flares with Stf/tf ~ 0.1 is 3 orders of magnitude less 
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than the energy that can be released during a smoother 
flare with Stf/tf ~ 1. The energy of the flare depends on 
E^°, the radial dependence of which has not been stud- 
ied for strongly magnetized ejecta. The key open question 
is how fast is the Poynting flux passed onto the forward 
shock. This has an important effect on the bulk Lorentz 
factor of the flow as a function of radius F(r) (see Zhang 
& Kobayashi 2005; Lyutikov 2006) and, therefore, on the 
(observer) timescales over which magnetic dissipation can 
result in powerful flares. A detailed study of the dynamics 
of the deceleration of ejecta with ao>l can settle these 
issues and yield valuable constraints on the model. 

Nevertheless, one can have a rough estimate of the 
energetics of a flare by assuming that the energy in the 
forward shock at some radius r is a large fraction of the 
energy that was initially carried by the ejecta E^°. In view 
of eq. 111(1 , it means that a fraction 

/<5e(^ f /t f ) 3 /a 2 (12) 

of the initial blast wave energy can power a flare of dura- 
tion Stf at (observer) time tf. As an arithmetic example, I 
set a 2 ~ 10 and e ~ 0.2, from which I find /<0.1(<Si f /i f ) 3 . 
One can go a step further and compare the energy of the 
flare with that of the prompt GRB emission by assum- 
ing that the latter is a fraction r\ ~ 0.1 of the energy 
of the blast wave. This means that a flare coming from 
late magnetic dissipation can reach a fraction ~ (5t[/t() 3 
of the prompt GRB emission. So the model may explain 
powerful and fast evolving flares at the same time. On the 
other hand, there is a clear prediction that fast evolving 
flares are less energetic than the smoother ones. 

2.4. Multiple flares 

Very often, the afterglow light curves are characterized 
by multiple flares. In the context of this model, this cor- 
responds to multiple reconnection regions. The physical 
scales of such region are constrained by eqs. © and 
and one can check that they correspond to a moderate 
fraction of the volume of a shell with 1/T opening angle 
that emits toward the observer. So, a single flare does not 
use up all the Poynting flux in the observer's light cone 
and repetition is possible. It is even likely that dissipation 
in one region can induce instabilities and magnetic dissi- 
pation in neighboring regions, which results in a sequence 
of flares. 

3. Conclusions 

The afterglow flares observed by Swift have been suggested 
to indicate late time activity of the central engine. As 
an alternative, I propose here that flares are produced in 
the deceleration phase of strongly magnetized ejecta. This 
comes out naturally from Poynting models where instabil- 
ities give partial magnetic dissipation (prompt emission 
and acceleration; see, e.g., Giannios & Spruit 2006). For 
high bulk Lorentz factor of the flow, the instabilities slow 



down because of the time dilation effect and the flow re- 
mains strongly magnetized in the afterglow phase. The 
deceleration of the flow due to interaction with the ex- 
ternal medium (and possibly the crossing of the reverse 
shock; Thompson 2005) revives the instabilities and leads 
to delayed magnetic dissipation. 

I have looked at the magnetic energy available in the 
flow to power such flares and how it depends on the length 
scales of the dissipation region. For a flare to have the ob- 
served rapid evolution, the scales of this region cannot 
be arbitrarily large, and this constrains the available en- 
ergy per flare. I estimate that flares caused by delayed 
reconnection can be powerful for both long and short 
GRBs, with isotropic equivalent energies up to a frac- 
tion - 0.1 • (Stf/tt) 3 of the blast wave energy. Although 
this estimate contains uncertainties related to the speed 
of magnetic reconnection and the deceleration profile of 
the ejecta [see eq. <|12f) ]. it indicates that powerful flares 
through delayed magnetic dissipation are possible and 
that smooth flares are expected to be more energetic than 
spiky ones. Multiple flares are expected in this model as a 
result of dissipation in multiple neighboring regions in the 
decelerating flow. Furthermore, inverse Compton scatter- 
ing of the flare photons that takes place in the forward 
shock can result in GeV-TeV flares as shown by Wang et 
al. (2006). 
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